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s a project tor an 
independent study- 
class, another student 
and I were implementing 
an IIR digital filter using a 68HC11 
microcontroller. We used MATLAB to 
design a second-order, low-pass IIR 
filter. MATLAB generated the transfer 
function that we implemented using 
the 68HC1 1. When we tested the fil- 
ter, we found the output was sporadic. 

Obviously, something was wrong 
with the filter, but we didn't know 
whether the problem was in the design 
or the implementation. To success- 
fully troubleshoot the filter, we first 
had to verify that the design was 
correct. If we could do this, then we 
knew our problem was in the imple- 
mentation, not the filter design. 

Because it was a student proiect, 
our method had to be inexpensive. We 
decided to use a spreadsheet to 
simulate and verify the filter design. 

IMPLEMENTING THE FILTER 
SIMULATION 

We used Microsoft Excel for the 
Macintosh to implement the simula- 
tion. Any standard spreadsheet can he 
used. However, to be most useful, it's 
best if the spreadsheet can plot graphs 
directly on the worksheet |see Figure 
1 1. To understand how the simulation 
works, let's first review digital filters. 

Digital filters sample an input 
signal and calculate the output value 



at specific, constant time intervals. 
The time between these intervals is 
the period. The sampling frequency of 
the filter is the reciprocal of the period. 
The characteristics of digital filters are 
based on the sampling frequency. 

Spreadsheets work well for digital- 
filter simulation because instead of 
being periodic in time, they're periodic 
in position. Each row in the spread- 
sheet can represent one sample of the 
input signal and the resulting calcu- 
lated output signal. 

PARTS OF THE SPREADSHEET 

The spreadsheet is composed of 
six parts: 

• sample column— serves as the 

timebase for the simulated filter 
and is used as a basis for the input 
and output plot. The input signal 
for the simulated filter is derived 
using the sample column. For most 
simulations, 100 data points are 
adequate. 

• input column — produces the simu- 

lated-input signal for the filter. The 
values in this column are computed 
based on the sample column and 
the values of the input signal 
frequency and sampling frequency. 
This is described in detail in the 
next section. 

• filtered column — contains the output 

of the simulated filter based on the 
input signal and the filter coeffi- 
cients. 

• signal and filter characteristics — 

specify the input signal frequency 
and the sampling frequency of the 
filter. They also specify the input- 
signal magnitude and any offset. 

• filter coefficients — come from the 

discrete-time transfer function used 
to describe the digital filter. These 
coefficients are used to calculate 
the values in the filtered column. 

• input and output plot — enables the 

input signal and filtered output 
signal to be viewed. The plot is 
produced by plotting the values in 
the input and filtered columns 
versus the sample column. 

PRODUCING THE INPUT SIGNAL 

The input signal is the most 
difficult part of the spreadsheet to 



Who says spread- 
sheets are just a 
financial tool? Steven 
proves otherwise. 
Cells easily 
accomodate input 
data and subsequent 
calculations of the 
output signal. Voila, 
your spreadsheet is 
a digital filter! 
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Figure ^-Microsoft Excel running on a Macintosh works well for filler simulation because it can display graphs of trie data on me same screen as trie spreadsheet itseil. 



implement. The input signal should be 
periodic and have characteristics of 
standard input signals (standard input 
signals include sine, square, and 
triangular waves). It also must 
demonstrate the correct relationship 
between the frequency of the input 
signal and the sampling frequency. 
Two basic calculations produce these 
results. 

To be periodic, the number of 
samples in one period of the input 
signal must be calculated. This value 
is determined by the ratio of the 
sampling frequency and the input 
frequency. 

Samples Per Period ■ ^ Frequency 
Input Frequency 

One cycle of the input signal must be 
completed in this number of samples. 

To correctly produce the input 
signal, the current position in the 
current period of the input signal must 
also be determined, as shown in Figure 
2. The position in the period |PP) is 
calculated by: 

PP-CurrentSarnp^l X^ZcT ) 

where % is the modulus operator. 
Eased on these calculations, the 
different input signals are produced 
(descriptions of some standard input 



signals follow|. Be sure to use absolute 
references for the input and sampling 
frequency, magnitude, and offset. Use 
a relative reference for the current 
sample. 

The sine wave is simulated using 
the S I N ( ) function found in standard 
spreadsheets. The signal is scaled 
based on the magnitude, rounded to 
the nearest integer, and then offset: 



Round j SIN ( P -i^ n g^f I I ■ 



The square wave is simulated 
using the I F ( ) operation and CO S ( ) 
function found in standard spread- 
sheets. When the output of the C0S( ) 
function is positive, the signal is the 
magnitude plus the offset. When the 
C0S( ) function is negative, the signal 
is the offset less the magnitude: 

\ Sampling Frequency J 
Magnitude + Offset 
ELSE 
Offset - Magnitude 

The triangle wave is simulated by 
two parametric equations. The signal 
starts at the positive magnitude and 
decreases to the negative magnitude by 
the middle of the period. In the second 
half of the period, the signal starts at 
the negative magnitude and increases 



to the positive magnitude. The I F ( ) 
operation determines which half of the 
period is currently being calculated 
(see Figure 3). 

The output of the filter is calcu- 
lated from the difference equation for 
the filter being simulated. The form of 
the difference equation varies depend- 
ing on the type of filter. When entering 
the formula to compute the output, 
you should use absolute references to 
the filter coefficients. Using this 
method, you have to write only one 
formula, which can be copied to all 
rows in the output column. 

USING THE SPREADSHEET 

This spreadsheet can be used to 
simulate IIR filters, FIR filters, and 




Figure 2— The current position in ttie current penod ot 
the input signal must oe determiner! tor the input signal 
to oe periodic. In this illustration, the value ot the input 
signal tor the fourth ot sixteen postlms is calculated 
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Genie 



Daia Genie offers a full line of test & measure- 
ment equipment that's innovative, reliable and 
very affordable. The 'Express Series' of stand- 
alone, non-PC based testers are the ultimate 
in portability when running from either battery 
or AC power. Data Genie products will be 
setting the standards for quality on the bench 
or in the field for years to come. 





HT-28 Express 



HT-14 **P re « 



DIGITAL IC TESTER ■ EPROM PROGRAMMER 



The HT-28 is a very convenient way 
oftesting Logic Cs and DRAM'S. Tests 
most TTL 74, CMOS 40/46 and DRAM'S 
4164-414000, 44164-441000. It can 
also identity unknown IC numbers on 
TTL 74 and CMOS 40/45 series with the 
"Auto-Search' feature. 
$189.95 



The HT-14 is one-to-one EPROM writer 
with a super fast programming speed 
that supports devices from 2732B to 
27080, with eight selectable pro- 
gramming algorithms and six pro- 
gramming power CVPP) selections. 
$289.95 




P-300 



PC INTERFACE CARD PROTECTOR 



The Data Genie P-300 is a useful device that allows you to quickly install add- 
on cards or to test prototype circuits for your PC externally. Without having to 
turn off your computer to install an add-on cards, the P-300 maintains com- 
plete protection for your motherboard via the built-in current limit fuses. 
$349.95 

MING 

Microsystems 

Division or MING ( 6 P. IHC^ 

17921 Rowland Street 
City of Industry, CA 9174B 
TEL : (818)912-7766 
FAX : (818) 912-9698 



Call tor a dealer near you. 

1-800-473-6606 



Data Genie products are backed by a full 
1 year limited factory warranty. 



CIPM MNG McotySemi. Al mswvW Daia Gente lege « eie^aetednc 



nt ol MNG Met 05ysf efm. 



analog filters that have been converted 
to discrete-time form. 1 suggest you 
make a template tor each type of filter 
and input, then begin experimenting. 
The description of three filters follow. 
Use these filter designs to test your 
templates. 

The following equation is the 
generalized transfer function for an IIR 
filter of order q, 

" i ^'" i X{2) = a fl + a 1 z- 1 + ... + a ^ 

In Figure 4, we see what convening 
this to a difference equation yields. 
The value of the difference equation is 
computed to produce the filter output. 
When calculating the difference 
equation, it's a good idea to use as 
many coefficients as the highest-order 
filter needs. When simulating lower- 
order filters, enter zeros for the higher- 
order coefficients. 

Note that the output is calculated 
from the current input and previous 
input and output values. For the filter 
to be causal, you need to use zero as 
the input for the samples prior to 
sample I, The number of these zero 
samples depends on the order of the 
filter you're simulating. Also, when 
you write the output formula, be sure 
to use relative references to the 
previous input and output values. 

After you've created the template 
for an UR filter, try simulating the 
filter shown m Figure 5a. This filter is 
a second-order, low-pass, Butterworth 
filter, designed using MATLAB for a 
sampling frequency of 1000 Hz and a 
cutoff frequency of 50 Hz. 

The following is the generalized 
transfer function for an FIR filter of 
order q. 

H|Z|«^^ = a fl + a 1 z- 1 +a 3 z- I +...+a (t z-< 

Converting this to a difference equa- 
tion yields: 

y|n)-aox|n) + a| x|n-l |+...+a^(n-q) 

Again, it's the difference equation 
that's calculated to produce the 
output. The difference equation is 
much simpler for a FIR filter, but FIR 
filters must be of a much higher order 
to have the desired characteristics. 
This makes FIR filters less practical to 
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Figure 3-Tfie lormuia 10 simulate a mangle-nave input signal uses ar, I - ) ooeralion to determine wnich Half ol lite 
currently being calculated. 



penod 



simulate using a spreadsheet, though it 
can be done. 

After you've created the template 
tor the FIR filter, try simulating the 
filter given in Figure 5b. This filter is a 
tenth-order notch filter designed with 



nant frequency of 100 Hz, and the 
output is the voltage across the resis- 
tor. The filter's transfer function is: 



Gfsl 



10 



10+ 0.2533s 



100000 



y|n 



b|>\ | n l-b,x | n - 1 1 +...+ b q x( 
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Figure 4— A difference equation is used to calculate the output ol an IIR met ol order q. 



MATL A3 for a sampling frequency of 
1000 Hz. it has a center frequency of 
125 FLz and a bandwidth ol approxi- 
mately 100 Hz. 

Analog filters can be simulated if 
they have been converted to a discrete- 
time form. The trapezoid rule is a good 
conversion to use to transform an 
analog filter from continuous-time 



Discretized using the trapezoid rule 
and a 1000-Hz sampling frequency, the 
discrete-time form is: 



H(Z) 



1-z- 2 



a) 
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Figure 5— Different kinds of filters— such as a second- 
order, low-pass IIR filter (a) and a tenth-order, notch 
FIR filler fly-can be simulated by inserting different 
coefficients. 

form to discrete-time form. The 
conversion is shown as: 

, 2|z-l) 
T|z+1) 

where T is the sampling period. After 
substituting for 5 in the continuous- 
time form and selecting a sampling 
period, simplify the expression and 
find the difference equation. Typically, 
the discretized transfer function has a 
difference equation like an IIR filter. 

As an example, try simulating the 
filter shown in Figure 6. It has a reso- 



'56.66- 91.3222-'+ 54.622 z- 2 

You can use the template you created 
for the IIR filter to verify the character- 
istics of this analog filter. 

EXPERIMENTATION 

After your templates work 
correctly, you can begin to experiment 
with the items listed below; 

• step response and settling time- 

make a template with a unit step 
input to view the step response and 
estimated settling time. 

• aliasing — watch for aliasing to occur 

when the sampling frequency is too 
low for the input signal. 

• output signal quality — note that the 

quality of the output signal de- 
grades as the input signal ap- 
proaches the sampling frequency. 




Figure h-One example of an analog filter that can be 
simulated is a band-pass filter with a resonant center 
frequency of 100 Hz. Trie output is across the resistor. 



• unstable filters— add an 
extra term to the denomi- 
nator of the example IIR 
filter. This will add an 
unstable pole to the filter, 
causing the output to 
"explode." (Fortunately, 
in the simulation, this 
iust means extremely- 
large output values, not 
physical damage.) 
• noise — use the RANDf 1 function in 
your spreadsheet to add simulated 
noise to the input signal, 

The only limits are the spreadsheet s 
capabilities and your own creativity. 

CONCLUSION 

Basic digital filters can be simu- 
lated using a spreadsheet, a tooi most 
people already have. Once you've 
created a set of templates, the simula- 
tions are easy to use, flexible, accurate, 
and best of all inexpensive. 

Getting back to my onginal 
problem, the spreadsheet simulation 
showed that the filter design was 
correct. The problem was in the 
hardware. After some debugging, the 
filter worked like the design. It was a 
beneficial problem, though. Now my 
colleague and I know that spreadsheets 
aren't just a financial tool, fg 

Steven Kubis is currently a senior 
technical writer with Gieat Plains 
Software in Fargo, North Dakota. He 
graduated in 1993 with a BSEE degree 
from North Dakota State University. 
He may be reached at skubis@ 
cogs.gps.com. 
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409 Not Useful 
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